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l/"") ' Abstract 

O ! 

, We present further properties of a previously proposed recursive scheme for parameteri- 

sation of n-by-n unitary matrices. We show that the factors in the recursive formula may be 
introduced in any desired order. The method is used to study the invariant phases of unitary 
matrices. The case of four-by-four unitary matrices is investigated in detail. We also address 
the question of how to construct symmetric unitary matrices using the recursive approach. 
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1 Introduction 

co . 

O ■ Unitary matrices play a central role in physics. For example, the Standard Model of particle 
Q I physics is defined by a SU(3) x SU(2) x U(l) symmetry group and many popular grand unified 
l/-) | models are again based on unitary symmetries. Indeed explicit representations of unitary matrices 
O ■ are often so badly needed that there is already a vast literature on the subject (see, for example 
^ ■ pQ and references therein). 

Recently we have presented a simple-looking recursive parameterisation of general n-by-n uni- 



^ ] tary matrices [2], applicable also, of course, to subcategories such as special unitary matrices and 
orthogonal matrices which are of great importance in physics. 

In recent publications Fujii and his collaborators [Sj have found that the parameterisation in 
looks interesting for constructing unitary gates for quantum computation but for that purpose 
^ . more study is needed. In this paper we present further results on the structure of the recursive 
parameterisation hoping that it will be useful for future applications. You use the method to study 
the "invariant phases" (to be defined below) of unitary matrices, by considering the symmetries of 
the recursive parameterisation. Subsequently, we give detailed attention to the case of four-by-four 
unitary matrices. 

It should be emphasised that all parameterisation of a general n-by-n unitary matrix are 
equivalent to each another. However, for a specific application a certain parameterisation may be 
more convenient than others. Therefore, it is important to provide new parameterisations, a topic 
which has been addressed by other authors as well. See, for example, [4 , a paper which contains 
an extended list of references and presents yet another representation of unitary matrices. 

2 The parameterisation 

A general n-by-n unitary matrix may be expressed as a product of three unitary matrices, 

lW=$W(a)y w $ (n) (^) (1) 
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where the matrices $ are diagonal unitary matrices, 



a) 



diag(e iai , e 



lOtl 



(2) 



$(/?) is defined analogously, the a's and /3's being real. We shall refer to $'s as external (pure 
phase) matrices. 

The matrix X^ n > has n 2 real parameters. The quantities a and (3 take care of 2n — 1 of these 
parameters because only the sums q.{ + j3j enter, where % and j run from 1 to n. The remaining 
(n — l) 2 real parameters reside in the non-trivial matrix which was the subject of the study 
presented in j2| and will be further investigated in this paper. For simplicity, whenever no confusion 
may arise, we refer to as the most general n-by-n unitary matrix leaving out the qualifying 
statement that this is only true modulus the external matrices $(df) and In [2 it was shown 

that the matrix may be written in the form 1 



n,n— 

where the A n>k are unitary matrices defined by 
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n.k 









In—k 



(3) 



(4) 
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Here I n - k is the unit matrix of order n — k. For k 
unitary matrix 

' 1 - (1 -c k )\A^ >< AW\ s k \AW > 
-s k <A^\ c k 



this unit matrix is absent. is a k-by-k 



(5) 



Here c k and s k stand for cosine and sine of an angle denoted by 9 k . Furthermore, \A^ > is a k — 1 
dimensional complex vector normalised to one, 



/ a[ k) \ 



\ A k) 



> 



\ 



M 
l k-i 



< a^\av >= l 



(6) 



/ 



and >< A^\)i 



a 



(*)(*)* 
a j 



. We shall refer to \A^ > as the characteristic vector of order k. 



The parameter counting was presented in j2] where it was shown that V^ n \ thus obtained, is 
the most general n-by-n unitary matrix, again modulus the external matrices $. The essential 
point is that \A k ^ > introduces 2{k — 2) real parameters and not 2(k — 1), the reason being that it 
is normalised and its overall phase can be absorbed into the definition of the external matrices <3>. 

To summarise, in this recursive parameterisation the n-by-n unitary matrix is represented by a 
product of n-1 unitary matrices, each with its own angle 9 and characteristic vector \A > while, for 
example, in the conventional approach in particle physics one would write the matrix as a product 

1 Note that the notations in this paper are a simplified version of that in [5] 
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of at least n(n — l)/2 matrices, these being Euler rotation matrices, (n-l)(n-2)/2 of them modified 
by phases (see, for example and references therein). Note also that in Ref. jl] the matrix 
is parameterised by a product of n diagonal unitary matrices interlaced with n — 1 orthogonal 
matrices. 



2.1 Reordering of the factors 

At the first sight, the recursive parameterisation appears highly ordered and rigid. In Eq.(jnj), the 
two-by-two structure is immediately followed by the three-by-three and so on. Actually, we may 
write these factors in any order we wish by observing that the order of factors in a given product, 
A n)T A n S) may be flipped as follows. For r < s we have 

An y rA n ^ s ^A n ^ r A n ^ s A^ nr }A n ^ r = A n s A n ^ r (7) 

where A' ns has the same form as A njS . The two characteristic vectors appearing in these matrices 
are related by a unitary rotation, 

\ A >(s) > = l(r)\ A (s) > (g) 

AM " ( T iL ) 

For the case s < r we have 

A nj rA ntS = An^^A^gAn^Ans) = A n:S A n ^ r (10) 

where now the two characteristic vectors are related by 

\ A n{r) >= AtW|AW > (11) 

Obviously, by inserting as many factors A\^A n j = 1 as needed in the recursion formula, Eq.(jHJ), one 
may move the factors around as one wishes. The upshot is that in the reordering process a factor 
of lower rank simply "tunnels" through that of a higher rank without being affected but induces 
a unitary rotation of the characteristic vector of the latter. Thus the ensuing parameterisation 
remains the most general one. Note that the angles 9k remain invariant under reordering. 

The recursive parameterisation looks highly asymmetric. However, using the above reordering 
procedure, one may construct manifestly symmetric unitary matrices (see Appendix). 

2.2 Further properties of the recursive parameterisation 

We wish to study, in more detail, the properties of the matrices A^ fc ^ in Eq.fjSj) as these are the 
building blocks of the recursive parameterisation. 

To begin with we simplify the notation, to avoid indices, and introduce a generic matrix A 
defined by 

k ^(l-{l-c)\A><A\ s\A>^ m 
— s <C yi c 
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Here, as usual, c and s stand for cosine and sine of an angle respectively. The angle itself will be 
denoted by 6. Defining 

Y=\AxA\ (13) 

we have that Y is hermitian and satisfies 

Y\A >= \A >, Y 2 = Y, trY = 1, detY = (14) 

where the vanishing of the determinant is, of course, only valid when Y is a matrix and not just 
a number as is the case when \A > is one dimensional. Following Fujii [Hj, we introduce a matrix 
G which generates A 

n = ( 

x i < A\ 

This matrix is hermitian and satisfies G 3 = G. A simple computation, using Eqs. (jl4j) . yields 

A = e ieG = 1 + isG- (1 -c)G 2 (16) 

This relation is reminiscent of the expansion of exponentials containing Pauli matrices a (in a 
short-hand notation, e t6a = c + isa). The essential point here is that the series expansion of e t9G , 
for arbitrarily G, terminates rapidly and does not continue for ever as the exponentials often tend 
to do. Furthermore we have 

trG = 0, trG 2 = 2, detA = 1 (17) 
Note also that, for a fixed G, the matrix A is Abelian with respect to 9, 

A(^)A(%) = A(0i + 6j) (18) 

and 

A- 1 (6)=A(-9) (19) 

To rewrite the recursion formula, Eq.(j2J, in terms of G and A we must attach appropriate indices 
to our generic G (or A) to distinguish the relevant factors. We introduce 

(( ^ A(k) >) \ 

&n,k= \t<A^\ J (20) 

V o 0/ 

where the required number of zeros have been added to make Q Hj k an n-by-n matrix. This yields 
that the factor A n ^ in the recursion formula Eq.(jnj) is given by 

A nM = e i6kG ">* = 1 + is k G n , k - (1 - c k ) Gl k (21) 
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3 Invariant phases of unitary matrices 



The invariant phases of a unitary n-by-n matrix V n are defined as those phases of the matrix that 
cannot be "removed" with any choice of the external phase matrices $ in Eq.fl^). These phases 
play an important role in particle physics as they are measurable quantities related to CP violation 
(for a review see, for example, Ref. [§]). 

Given a unitary matrix, the simplest way to detect the presence of invariant phases in it is to 
construct 

(a(3;jk) = Im(V aj V pk V: k V^) (22) 

where we have suppressed the superscript n. The symbols a, (3 and j, k now refer to rows and 
columns of the matrix and the indices are not summed. These imaginary parts are manifestly 
invariant under multiplication by the external phase matrices. Therefore if any of them is nonzero 
that would be a signal of the presence of a nonremovable phase in the matrix. We refer to these 
imaginary parts as invariant phases of the matrix instead of calling them invariants of the matrix 
that contain nonremovable phases. One may easily construct higher order invariants, containing 
properly chosen, six or more, elements of the matrix but these are in general reducible to the 
above set unless the matrix would have vanishing elements. For example, for V@j ^ 

W/sWW = (23) 

= ^{(aPJk)(J3 7 Jl) + (aPJk){P'r 1 jl)} (24) 
\Vpj\ 2 

where none of the indices is summed and 

(a/3; jk) ee Re[V aj V pk V* k Vl\ (25) 

These real parts are also invariant under the action of the external matrices. The above reduction 
would not work if V a j = but then the analysis is much simpler to begin with (see below) as the 
matrix contains fewer invariant phases. 

Returning to the simplest invariant phases, there are altogether [n(n— 1)/2] 2 quantities (a/3; jk), 
because these are antisymmetric under the interchange of the row indices, a <-> (3, as well as under 
the interchange of the column indices, j <-> k. However, we know that the most general has 
"only" (n — l)(n— 2))/2 independent invariant phases. One may therefore look for (n — l)(n— 2))/2 
independent (a(3; j7c)'s and use them as a basis for expressing the remaining ones. 

As mentioned above, the invariants in Eq. (j22|) play an essential role in the n-family version of 
the Standard Model of particle physics as they are measurable quantities related to CP-violation. 
For the case of n = 3 there is only one such quantity 

(a(3;jk) = J )^ e ja pe ijk (26) 

The row and column unitarity conditions for a three-by-three unitary matrix define six triangles. 
One may show that [7] all these triangles have the same area and this unique area equals J/2. 
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For n = 4 there are 36 possible invariants (a(3;jk) but only three independent ones. In [Hj 
an attempt was made to find an appropriate basis and carry through the above programme. The 
treatment of this issue is much simpler in the recursive parameterisation, as will be shown in the 
next section. 



3.1 Invariant phases of four- by-four unitary matrices 

For n = 4 we have from Eq.(J3J) 

V® = A 4t2 A 4>3 A 4>4 (27) 
where each factor comes with its own 9 and characteristic vector \A>. We denote the latter by 



2) 



\A^ >-- 



X2 



\A& >-- 




(28) 



remembering that x's and y's are complex numbers and < A^ k '\A^' >= 1, k = 2,3,4. We shall 
now spell out this four-by-four matrix in order to exhibit its symmetries in a manifest fashion. 
This will also enable us to understand the general case of n-by-n matrices. Eq.Q yields 



I c 2 s 2 \ / 1 - (1 - c 3 )xix\ -{l-c^xxx^ s 3 xi 0\ 



X 



s 2 c 2 
10 
\ 1 / 

/ 1 - (1 — 04)2/12/* 
-(1 - 04)2/22/1 

-(1 - C 4 )2/32/t 

V -S42/* 



'1 - c 3 )a;2a;* 




1 - (1 - C 3 )X 2 X2 S3X2 

S3X2 c 3 





-(1 - C 4 )2/l2/2 -(! - C4)2/l2/3 S 42/l \ 

1 - (1 - C4)2/22/2 -(l-C4)2/22/3 S 4 ?/2 

-(1 - 04)2/32/2 1 - (1 - 04)2/32/3 S 42/3 

-S42/2 - S 4?/3 C 4 / 



1/ 



(29) 



We now focus on the symmetries of this matrix. By symmetries we mean transformations that 
leave invariant modulus the external matrices $ in Eq.([T|). A simple inspection shows that 
this matrix has two such symmetries, denoted by Si and S* 2 and defined by 



Si 



So 




Mf>Z 




2/3 



e-^ 2 y 3 



(30) 
(31) 



where 02 and 03 are arbitrary phases. The indices are to remind us of the dimension of the 
corresponding vector. Therefore the three independent phases in can be chosen to be: 



ui = <j>{x 2 ) - 4>(xi) 
^2 = 0(2/2) - 0(2/i) 
^3 = 0(^2) + 0(2/3) 



0(2/2 



(32) 
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Here <f>(xj) and <f>(yk) denote the phases of the corresponding parameters. These phases are not 
invariant under the above symmetries and thus can't appear as independent entities in computa- 
tions of invariants of . Allowed to appear are the u;'s or any combination of them because 
these are invariant under the action of both Si and S2. We may, if we so wish, use the symmetry 
Si to rotate the phase of x 2 to zero, and then employ the symmetry S2 to do the same with the 
ensuing y 3 whereby x 2 and y 3 may be taken to be real and say positive. The invariant phases in 
this "frame" are then (f>(xi), 0(yi) and 0(1/2)- These constitute the maximum number of indepen- 
dent phases that can possesses. By imposing further relations on the angles or the x's and 
y's this number could be smaller as shall be considered further below. 



3.2 Generalisation to larger n 

Going one order higher to n = 5, we have to introduce the relevant characteristic vector 

\A^ >= 



Z 2 
Z 3 



The corresponding matrix has three symmetries given by 
Si : 



S 2 



S 3 




— 100 

e V2 y 3 , z 3 



—i<p3 



In this case there are six invariant phases. These may be chosen as 



UJi 


= <f>(x 2 ) ~ <f>(xi) 




UJ 2 


= 0(1/2) - 0(yi) 




L0 3 


= <P{z 2 ) - (f)(zi) 




U4 


= 0(X 2 ) + 0(l/3) - 


0(l/2) 


L0 5 


= 0(2:2) + 0(^) - 


0(^2) 


L0 6 


= 0(l/3) + 0(^) - 


0(^3) 



(33) 



e- l<p2 z 3 



(34) 
(35) 

(36) 



(37) 

As before, we may use Si to remove the phase of x 2 followed by S 2 and S3 to rotate the ensuing 
1/3 and z 4 to be real and say positive. The invariant phases are then the phases of Xi, yi, y 2 , 
zi, z 2 and z 3 . Note that, from the very beginning we chose the angle 9 2 not to be accompanied 
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by a phase, i.e., \A^ >= 1. One can, of course, leave the phase of \A^ > arbitrary. This will 
introduce an extra symmetry which we have not bothered to write down as it is trivial. 

The above procedure may be generalised to arbitrary order n. Without loss of generality, we 
may take, for example, the last component of all the characteristic vectors \A^ kS} > to be real. The 
invariant phases are then the phases of the remaining components. For an n-by-n matrix there 
are then 1 + 2 + .. + {n — 2) = (n — l)(n — 2)/2 such independent phases as expected. 



3.3 The "panel" approach to invariant phases 

Another approach to constructing the invariant phases of a matrix is to consider the latter as a 
lattice. For the case of n — 4, considered from now on, the matrix can be visualised as shown 



Pll 


P12 


Pl3 




• 


• 


P21 


P22 


P23 




• 


• 


P31 


P32 


P33 



(3* 



The bullets denote the sites where the matrix elements are situated. For example the bullets on 
the first row stand for Vn, V12, V13 and Vu and so on. The P's denote minipanels of the matrix, 
to be described here below. The invariants of interest to us are again 

(a(3;jk) = Im[V aj Vp k V* k V^] (39) 
(aftjk) = Re[V aj V pk V: k V^} (40) 

As mentioned before, there are 36 quantities (a/3;jk) (six possible combinations of a, (3 multiplied 
by as many combinations of j, k) and we are looking for a set of three of them such that all the 
others can be expressed as functions of them and the real parts (a/3; jk). This problem was treated 
long time ago : 8j and was found to be rather involved. Here we provide some simplification. Nine 
of these invariants (the nearest neighbours) are explicitly exhibited on our lattice. Their analytic 
form may easily be read off from their locations. For example, 

Pll = VllT^22^12^21 
Pl2 = ^12^23 ^13^22 
-^22 = ^22^33^23^2 
-P32 = V32 V^K^V^ 

and so on. Furthermore 

Pab = Rab+iJab (41) 

where R and J denote the real and imaginary parts of the corresponding P. Thus the imaginary 
parts, in the notation employed earlier, are given by 

Ju = (12;12), J 12 = (12;23) J 22 = (23 ; 23) J32 = (34; 23) (42) 
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and so forth. Similar expressions may be written down for the real parts. 

Suppose that none of the matrix elements vanishes. Using unitarity conditions we then find 



•A3 - (1 + . 11 l2 )Jn - |T . \ 2 r n -h\ 



■^"^33 ~^~^23 
^13 - (1+ , T/ , 2 )^23 = , T/ 77 , 2 ^33 



^31 - (1 + , T/ 7T 12)^21 - , T/ ^ , 2 ^11 
|V2lV22r I V2I V22 I ^ 

■^"^33 "^^32 

|V3 3 V4 3 r 1 V33V43I 

-^^22 "^^32 

•^12 - , T/ T/ ,2^32 = (1 + , T/ 77 ,2)^22 
V32V33M K32K3." U 



/ 32> / 33| I ^32 V33 I 

^21 - 777^^23 = (1 + 777^)^22 (43) 
I ^23 ^33 1 



In principle, we may take the set Jn, J22, J33 to constitute a basis and determine the remaining six 
J's in terms of them. These equations are rather complicated and need further thought concerning 
special cases. For example, if the matrix is symmetric we only have three equations but, of course, 
also only three unknown, say J12, J13 and J23. We are not allowed to divide by vanishing matrix 
elements and so forth. Here below, we shall consider a simple and yet nontrivial example to 
demonstrate the technique and to compare it with the recursive approach which is much simpler 
and does not require thinking about the possible pitfalls. 



3.4 A simple example 

For simplicity we consider the case where two of the elements of the four-by-four matrix are zero, 
and where these elements are neither on the same row nor on the same column. All other elements 
of the matrix are assumed to be nonzero. Without loss of generality we may take the two vanishing 
elements to be V14 and V41 which in particle physics would correspond to the case where the mixing 
of the first and the fourth families is negligible. Our lattice, with its nine minipanels, now looks 
as follows 

• • • © 



J 


J 





• 




• 


J 


J + J' 


J' 


• 




• 





J' 


J' 


• 




• 



Here the 0's indicate where the vanishing matrix elements are situated and we have defined 

J= J n = (12,12), J'= J 33 = (34,34) (45) 
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A simple computation, using unitarity relations, gives the imaginary parts of the minipanels as 
marked in the lattice. Computing all the imaginary parts, we find that 19 of the 36 invariants 
(a/3,jk) vanish. The nonvanishing ones, in addition to J and J' defined in Eq. (145(1 . are 

- (12,13) = (12,23) = -(13, 12) = (13, 13) = -(13,23) = (23, 12) = -(23,13) = J 
-(23, 24) = (23, 34) = -(24, 23) = (24, 24) = -(24, 34) = (34, 23) = -(34, 24) = J' (46) 

and 

(23, 23) = J + J' (47) 
as exhibited in the corresponding panel. Moreover, we find 

J' 2 i V24V34 1 2 _ I ^42^43 1 2 



J 2 , v 21 v 31 l >vM 

, I V24 1 2 + I V34 1 2 2 _ I V42 1 2 + I V43 1 2 . 2 



(4S 



It is amusing to note that the unitarity relations for the above matrix define eight triangles. Using 
the method in [7j one finds that four of these have each an area equal to J/2 while the area of the 
other four is J'/2. 

We would now like to compute J and J'. For this purpose we turn to the recursive parame- 
terisation. It turns out that the calculations are simpler is we take V34 = V43 = instead of the 
above choice Vu = V41 = 0. The two choices are equivalent as they are related to one another 
by interchanges in rows and columns. This amounts to a relabeling of the matrix elements which 
obviously can't affect the results. After finishing the computations we can simply revert to the 
former case by interchanging rows one and three as well as columns one and three. 

In the recursive parameterisation, Eq. (j35|l . the conditions V34 = V43 = give 

2/3 = xij/J + x 2 y* = (49) 

These condition tell us that l^l = \x 2 \, I2/2I = Nil an d that X\ and yi are relatively real in the 
frame where x 2 and y 2 are taken to be real. Therefore, there is only one invariant phase, in this 
example. We introduce the lattice again 



J 


J 


• 




J 


J 


• 










• 





(50) 









0's indicate where the vanishing matrix elements are situated. Furthermore, we have exhibited 
the imaginary parts of the minipanels starting with the definition 

J = (12, 12) (51) 
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Taking into account the permutations, we find that J here is identical with our previous J that 
we wanted to compute. Using the recursive parameterisation we find 



J = C2C3C4S 2 sllm(XiX 2 ) 

J' = -c 2 C' i c i S2s\lm(x\x2) (52) 

Thus 



.r _ s j 

J s 2 



(53) 
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For comparison, note that for a general three-by-three matrix (which we may obtain from Eq. 
by putting 6 4 = 0) the unique invariant is given by j( 3_ / am ) = c2C 3 S2s\Im(x\x2) ■ 

The recursive parameterisation allows us to compute all the imaginary parts for the most 
general case, i.e., irrespectively of whether the matrix has zeros or not. We find, for the general 
four-by-four matrix, parameterised as in Eq. ()29j) 

(34,34) = c 3 CiS 3 sl\y 3 \ {\x 2 y 2 \sinuj 3 + Ixxyxlsinius + u 2 - uJi)} 

(34,24) = c±s 3 s\\x2y 2 \{s 3 \x 1 y 1 \sin{u 1 - u) 2 ) - c 3 \y 3 \simo 3 } (54) 



where the angles u>j are as defined in Eq. (j32|) . We thus see that, as expected, only the invariant 
phases appear in these relations. We do not quote the remaining imaginary parts (af3,jk). The 
important point is that all of them are functions of \xj\, \yj\ and the three u/s, as expected. 



4 Conclusions 

In this paper, we have presented further properties of the recursive parameterisation of unitary 
matrices proposed in |2j, where the matrix is written as a product of n — 1 matrices each with 
its own angle 9 and characteristic vector \A >. We have found that the factors in the recursive 
formula may be introduced in any desired order. 

Encouraged by the convenience of the recursive method, we have taken a fresh look at the issue 
of invariant phases of unitary matrices. After having exhibited the symmetries of the parameteri- 
sation, we have shown how the invariant phases of n-by-n matrices can be identified. Subsequently, 
we have paid particular attention to the case n = 4 and have compared the results with those of 
an earlier approach based on "panels" of the matrix. 

The recursive parameterisation has some really nice features because in some cases it allows 
the "new physics" to be introduced in a gentle manner through the last factor in the recursion 
formula, a topic which we are currently studying. 

In an earlier study [Oj, we found that there is a parameterisation that allows one to introduce, in 
a simple way, any desired angle of any of the so called unitarity triangles as one of the parameters 
in the quark mixing matrix for three families in the Standard Model of particle physics. The same 
parameterisation allowed us to choose the expansion parameter in this matrix to be A 2 instead 
of A that one usually uses [5]. Indeed A is not so small (A = 0.2). Therefore, the recursive 
parameterisation may be convenient whenever expansion in the above parameter is required, for 
example in model building or for construction of quark and lepton mass matrices. It turned out 
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that the parameterisation found in [Hj, with the above nice features, is indeed nothing but the 
order n = 3 version of the recursive parameterisation discussed in this paper and in 

Finally, in the appendix of this paper, we deal with the question of how to construct manifestly 
symmetric unitary matrices in the recursive framework. 



5 Appendix: Symmetric unitary matrices 

We write the symmetric unitary matrix in the form 

X (n)sym = $ (n)^y( B )i|im $ (n)^ (55) 

requiring V^ sym to be symmetric, as indicated by its superscript, and that the external matrices 
be the same (see Eq.flQ). A general symmetric unitary matrix has n(n + l)/2 real parameters. 
The external matrix <3> takes care of n of them. Thus n(n— l)/2 real parameters reside in V^ n ' 8ym . 

For a general V^ n \ we have (see Section F2.2j) that the factors in the recursion formula Eq.fjHJ) 
may be written as 

A nM = e i9kG "> k (56) 

where the generating matrix G n ^ is hermitian. In order to obtain a symmetric A n ^ we must 
impose the additional requirement that the generating matrix be symmetric. This means that the 
corresponding characteristic vector is purely imaginary. For example, for k = 2 we obtain 



y(2)sym = j^m = ^ Q (57) 




A sym 



(5f 



and for k = 3 

1 — (1 — Cz)x\ —(1 — C%)XiX2 IS3X1 

— (1 — c^)x\x2i 1 — (1 — 03)^2 IS3X2 
253X1 IS3X2 C3 

\ J n _ 3 / 

We have put \A >= i\x >, where the x's are real. Thus, the construction of symmetric factors 
in the recursion formula is a trivial task. But, of course, the product of these factors will not 
be symmetric. This defect is easily remedied by invoking the reordering procedure described in 
Section 12.11 in which we showed that the reordering of the factors in the recursion formula only 
amounts to a redefinition of the characteristic vectors. Therefore, we may write 

T/(n)sj/m _ Asym Asym A sym Asym /.sym Asym Asym /r Q \ 
v — ^n,2 ^,3 ••• /1 n,n-l /1 n,n /± n,n-l--- /i n,3 -^1,2 l dy J 

y(n)sym obtained is manifestly unitary and symmetric. We must now count the number of its 
independent parameters. Each order k introduces k — 1 real parameters, these being the angle 9k 
and k — 2 components of the corresponding characteristic vector (one component being redundant 
because the vector is normalised). Therefore the total number of parameters in V^ nS>sym is 

n 

J2(k-l) = n(n-l)/2 (60) 

k=2 
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as expected. Adding into this number the n parameters coming from the external matrices amounts 
to the total of n(n + l)/2 real parameters, as required. 

Note that it would be somewhat more elegant to call the angles in the above factors 6^/2 
instead of 9k, except the angle 9 n of the factor A sy ™ . The reason being that A sy ™ appears only 
once while the others appear twice. 

The chain in Eq. (|59j) looks long and perhaps a bit frightening. However, if needed for practical 
applications, it can be somewhat simplified as we shall now describe. 

Consider the case n = 3, where we introduce 



y(3)sym = y (2)sym (Q 2 /2)A^ Sym V^ Sym (9 2 /2) 

where V^ sym is as defined in Eq. (157)1 . for n = 3. Multiplying the factors, we find 



(61) 



c 2 — (1 — c 3 )u\ is 2 - 
r i:iN ""' = I is 2 -{l-c 3 )u 1 u 2 c 2 
is 3 u\ 



(1 - c 3 )miM 2 is 3 ui 
(1 - c 3 )u 2 2 is 3 u 2 
is 3 u 2 c 3 



(62) 



Here 



u\ = c 2 x\ + is' 2 x 2 , u 2 = c 2 x 2 + is' 2 x\ (63) 

cos(9 2 /2) and s' 2 = sin{9 2 /2). Moreover, the x's are as introduced in Eq.fjHHJ)- 
The essential point is that we may put aside the question of the origin of the u's and their 
relationship with the x's and simply consider them as our new variables, as two complex numbers 
that satisfy 

\ui\ 2 + \u 2 \ 2 = 1 (64) 
Going to the next order, n = 4, we may use the identity 



c' 2 x 2 + is' 2 xi 



4 



A sym a sym a sym a sym a sym 
^2 ^3 ^4 ^3 A 2 



A' 4 sym 



j^sym j^sym j^syrn Uj^ s V™-\ 



-lj^sym^jyym^ 



-lj j^sym j^sym j^sym 



{A sy m) 



sym\—l 



(65) 



Here V^ sym is as found in Eq.()62|). From our earlier results, we have 



jyym 



( 1 - (1 - c,)y\ 
-(1 - c 4 )yiy 2 
-(1 - c A )y x y 3 

\ ism 



-(1 - c 4 )ym 
1 - (1 -c 4 )2/| 
-(1 - c A )y 2 y 3 

ism 



-(1 - c 4 )yiy 3 ism \ 

-(l-c 4 )y 2 y 3 is 4 y 2 

1 - (1 - c A )yl ism 
ism c 4 J 



(66) 



where y's are real. Therefore, we may immediately write down the factor A'^ ym , without having to 
do any calculations. The first two components of the vector y get "rotated" but y 3 is untouched. 
We find 



psyml 



( C 2 - (1 - c A )v\ 
—is 2 - (1 - c 4 )viv 2 

-(1 - Ci)viV 3 

ism 



-is 2 - (1 - c 4 )vxv 2 
c 2 - (1 - c 4 )v 2 
-(1 - c A )v 2 v 3 
ism 



-(l - c 4 )vm ism \ 

-(l-c A )v 2 v 3 ism 

1 — (1 — c A )v 3 ism 
ism c 4 / 



(67) 
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where 

vi = c 2 yi - is' 2 y 2 , v 2 = c' 2 y 2 - is' 2 y x , v 3 = y 3 (68) 

Again the vector v has unit norm and we may, as before, forget about the y's and just use i>'s, 
keeping in mind that V\ and v 2 are complex numbers. 

As a final example, we wish to compute the quantity J for the case of a three-by-three sym- 
metric matrix, where 

(a(3;jk) = J ^ e -ra^ijk (69) 

and V^ sym is as given in Eq.(|62j). A glance at this matrix yields 

J = c 2 c 3 sllm(u 2 ) 2 = C2C3S2S3X1X2 (70) 

This resembles our earlier result in Section I3~H where we found J^f am ) = c 2 c 3 s 2 s\Im{x\x 2 ) . The 
meaning of the x's in the two cases are, of course, different. Note that a general has four 
parameters while V^ sym has one less. Equation (|70p is telling us that the symmetry requirement 
does not remove the invariant phase of the matrix. In the language of Euler rotations, where 
would be parameterised with three rotation angles and one phase, the requirement that the matrix 
be symmetric keeps the phase but removes one of the rotation angles. Note that the new phase 
and angles will be functions of the former phase and angles. 

References 

[1] R. Gilmore, "Lie groups, Lie algebras, and some of their applications" (Wiley- Interscience, 
1974); 

H. Georgi, "Lie algebras in particle physics: from isospin to unified theories" (Reading, Mass., 
1982) 



[2] 
[3] 

[4] 
[5] 
[6] 
[7] 
[8] 
[9] 



C. Jarlskog, math-ph/0504049, to be published in Jour. Math. Phys. 



K. Fujii,|math-ph/0505047; 



K. Fujii, K. Funahashi and T. Kobayashi quant-ph/0508006 

P. Dita, J. Phys. A: Math. Gen. 36 (2003) 2781; |math-ph/0103005| 

Review of Particle Physics, S. Eidelman et al, Phys. Lett B592 (2004) 1 

C. Jarlskog in CP Violation, Ed. C. Jarlskog (World Scientific, 1989) p. 3 

C. Jarlskog and R. Stora, Phys. Lett. B208 (1988) 268 

C. Jarlskog, Phys. Rev. D36 (1987) 2128 

C. Jarlskog, Phys.Lett. B615 (2005) 207-212 Qhep-ph/0503"T99"l ); 
see also C. Jarlskog, hep-ph/0504012 



14 



